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I . INTRODUCTION 



Today; naval architects are greatly interested in the design of 
super tankers. However; many difficulties arise in design when the 
strength of the super tanker is considered from an economic point of 
view. The transverse strength problem is apparently one of these diffi- 
culties . 

There are several ways in analysing the transverse strength prob- 
lem. In each method; the two-dimensional approach is simpler; but the 
three-dimensional approach is so complicated that hardly no one wishes 
to use it. Therefore; some approximate is devised and sufficient results 
are achieved. But still some architects wish to know the three- 
dimensional effect on transverse strength analysis. 

Some prominent naval architects suggested solving the three- 
dimensional strength calculation but most of them did not apply numerical 
computations to existing complicated ships ; for even though the computer 
is available; it is still no easy task. So many variables exist it is 
very tedious work to put them into computers . 

In this paper we are going to consider these tedious problems by 
using FORTRAN 7090* 

As previously mentioned; there are many ways to solve the three- 
dimensional transverse strength problem and each one has its own merits; 
although they will give almost the same level of accuracy. 

Since some approaches are so complicated and give the same level 
of accuracy as more simple approaches; they become useless. Therefore; 
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we are going to deal with a well-known and more simple approach- -the 
slope-deflection raethod--to attack these three-dimensional calculations 
and examine its effect. 

In large tankers especially^ each span has such deep depth- 
breadth ratios that we suppose shear force contributes more than usual^ 
therefore^ we take into account shear force effect in slope deflection 
equation . 

Eventually^ we will check our computed results by means of 
slope-deflection method with one of the fixed point methods which was 
given in our class by Professor Schade. This fixed point method differs 
in assumption from the slope-deflection method^ but we can make this 
fixed point assumption comparable with our slope-deflection assumption 
by restricting our loading conditions. Although we are dealing with the 
fixed point method^ it will not give us shear force effect and three- 
dimensional effect. We can see^ therefore^ only the differences between 
those two: the slope-deflection method with shear force effect and the 

fixed method without. They will lead us to some interesting conclusions. 



II. BASIC SLOPE-DEFLECTION EQUATION CONSIDERED WITH SHEAR EFFECT 



(l) Simple beam 
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Energy equation of the beam 
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Applying Castigliano theorem 
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(2) Simple beam with bracket at both ends 
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Rb = % - 



\ = '^C ■ "c - Va 

''b = «D - " «b(^B 
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Where 




also 



®C “ ®A ©D “ 

/o = i C4i - - W 



i = A, B 



Using all of these results and the result of (l), we can easily get 
the following basic formula using the following notation. 
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Where 
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and the constant terms are: 
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^mn >“ sum of the load on M'N' 

W'mn > sum of the load on MM' 

W' ^ ^ sum of the load on N'N 

nra 

Smn ’ s'mn ■’ sJ'nm distance from left end edge to 

the center of load for each span. 



For the purpose of our calculations, we only need two kinds of 



loading spans. They are shown on Table I and Table II. 
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Table 1 
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Table 2 
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III. SETTING UP SIMULTANEOUS EQUATION 

A. Assumption 

Once the designer has made the decision to consider the trans- 
verse frame of a ship^ he is then faced with the task of deciding on 
simplifying assumptions. The main objective of these assumptions is to 
create a fictitious structure and condition of loading^ which will re- 
duce the amount of labor in performing the structural analysis without 
compromising too severely the validity of the numerical answers. 

(1) One of the factors a designer must consider is the degree of 
dimensional similitude he can retain in the structure he selects for his 
computational work. For example^ he must decide how much of the shell 
plating should be included in the analysis of a transverse frame as ef- 
fective material. We are going to use Effective Breadth Conception in 
Ship Structure Design ^ published in 195^ tiy Professor Schade. 

There are several other assumptions to be applied in this paper. 

(2) The relative angle initially formed by the tangents to the 
elastic curve at a point where two or more members are joined are 
assumed to remain unchanged after loading; that is^ if two members are 
originally joined at right angles^ the angle formed by the tangents to 
the elastic curve after the loads are applied will remain ^ 0 ^ . 

(3) Joints of 1-1 (or 9-2) and 4-2 (or 13-2) are fixed in space 
but free to rotate in three-dimensional analysis and additionally 9“1 
(or l4-2)^ 15“2 (or 13 -l)^ l4-l and 15-1 are fixed in space in two- 
dimensional analysis (see Fig. 2-2). 

B. Simultaneous Equation 

Before setting up simultaneous equations we must investigate the 
number of unknown variables. From the basic slope equation we can take 
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up the following variables as unknown variables . 

1. 91(1) = 02(9) = X(l) 

2. 01(9) = 01(5) = 02(14) = X(2) 

3- 02(1) = 02(10) = 91 ( 2 ) = X( 3 ) 

4. 02(5) = 91(10) = 01(6) = X(4) 

5. 02(2) = 02(11) = 01(3) = x( 5 ) 

6. 92(6) = 01(11) = 01(7) = X(6) 

7. 02(3) = 02(12) = 01(4) = x( 7 ) 

8. 92 ( 7 ) = 01(12) = 01(8) = x(8) 

9 . 92(4) = 02(13) = x( 9 ) 

10. 02(8) = 01(13) = 02(15) = X(10) 

11. jz((l) = s?i(5) = x(ll) 

12. jzi(2) = sz((6) = X(12) 

13. = sz^(7) - X(13) 

14. </>{h) = sz((8) = X(i4) 

15. H 9 ) = iz^(lO) = jzi(ll) = sz(( 12 ) = szi(i 3 ) = X( 15 ) 

16. i;i(l4) = X(l6) 

17. ^{ 15 ) = X(17) 

Now we have I7 unknown variables in one transverse frame. 
Therefore we have to have at least I7 simultaneous equations for each 
frame . 

1. Ml(l) + M2(9) = 0 

2. Ml(9) + Ml(5) + M2 (i4) = 0 

3. M2(l) + M2(10) + Ml(2) = 0 

4. M2(5) + Ml(lO) + M1(6) = 0 

5. M2(2) + M2(ll) + m.(3) = 0 

M2(6) + m.(ii) + Mi(7) = 0 



6 
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7. M2(3) + M2(12) + - 0 

8. M2(7) + M1(12) + KL(8) = 0 

M2(4) + M2(13) = 0 

10. M2(8) + ML(13) + M2(15) = 0 

11. 01(6) + 01(2) - 02(1) - 02(5) = 0 

12. 01(7) + 01(3) - 02(2) - 02(6) = 0 

13. 01(8) + 01(4) - 02(7) - 02(3) = 0 

14. X(ll) X SPM(l) + X(12) X SPM(2) + X(l3) x SPAN(3) + X(l4) 

X SPAN(4) = 0 

For the three dimensional analysis we can set up the following 
three more equations . 

15 . Yg = X( 15 ) X SPM( 13 ) + X( 17 ) X SPAN(15) 

16 . Ip = X( 15 ) X SPM(9) + X(l6) X SPM(14) 

17 . = X(15) X SPAN(13) 

Here and are deflection of bottom and deck longitudi- 

nals and longitudinal BHD. We suppose that longitudinals are only acted 
by the concentrated load at joints with transverse frames . Then acting 
load and are expressed as follows: 

Pg = 201 ( 15 ) 

Pd = 201(14) 

P^ = 01(9) + 01(10) + 01(11) + 01(12) + 01 ( 13 ) - 02 ( 15 ) - 02(14) 
For the two dimensional analysis we need only to solve l4 simultaneous 
equations with unknown variables l4 by setting simply X(l5)> X(l6) and 
X(17) to zero in the above equations. 

In this paper ^ for the practical numerical computation^ we will 
take an existing Japanese built Tanker as an example. She has 47^000 
tons Dead Weight and one of the tankas length is 503-936^^' and height 
604 . 329 ”^ and its breadth is 1188.97^^^’* of the web frame 
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scantling is shown in Fig. 1 and three weh frames placed in a length of 
tank. 







'SDTTQtl.....i>h^ J3HQXt^eSS__j^33£_ 



FIG. 1 TRMS. WEB FRAME 
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IV. PRELIMINARY CALCULATION 

A. Section modulus calculation 

From the effective breadth conception curve CASE (lll)^ ve can 
easily find out p : Effectiveness ratio. From these data we can calcu- 
late moment of inertia and neutral axis. 

The results are as follows: 



SPM No. 


F 


MOMENT OF INERTIA 
OF THE SECTION 


RATIO OF 
INERTIA 


DISTANCE FROM N.A. 
TO FREE FLANGE 


1 


0.570 


30828.851 


0.913 


47.240 


2 


0.350 


25839.827 


0.794 


42.339 


3 


0.350 


25839.827 


0.794 


42.339 


4 


0.440 


45197.227 


0.822 


41.335 


5 


0.635 


23604.205 


0.891 


42.189 


6 


0.410 


20716.733 


0.780 


38.867 


7 


0.410 


21014.376 


0.784 


39.184 


8 


0.485 


60847.290 


0 . 8 l 4 


33.413 


9 


0.870 


22181.710 


0.985 


44.287 


10 


1.000 


3083.156 


1.000 


16.259 


11 


1.000 


3656.705 


1.000 


16.377 


12 


1.000 


5282 . 347 


1.000 


16.748 


13 


0.870 


84741 . 346 


0.974 


55.841 


14 


0.735 


24303.928 


0.963 


43.268 


15 


0.735 


60303.312 


0.945 


57.421 



Table 3 




I 



Ik. 



In inches 



SPAN 


SPAN LENGTH 


e 


_ ^1 


No. 


L 


^2 




198.454 


122 . 859 


58.367 


1 


17.228 


o 


125 .984 


90.967 


17.372 


c. 


17.645 




125.984 


89.928 


17-645 


3 


18.411 


4 


137.987 


40.587 


17-843 


79-554 


5 


198.454 


120.428 


60.574 


17.452 


6 


125.984 


90.536 


17.596 


17.892 


7 


125.984 


89 . 456 


17-857 


18.671 


8 


137-984 


39.287 


17.955 


80.745 




301.006 


176.205 


61.037 


9 


63.764 


10 


301.006 


196.093 


52.051 


52.862 


11 


301.006 


196.840 


51.670 


52.496 


12 


301.006 


198.488 


50.850 


51.668 


13 


301.006 


190.620 


54.579 


55.807 


14 


280.539 


233.377 


11.062 


36.100 


15 


280.539 


236.469 


12.451 


31.619 



Table 4 
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Fig. 2-1 Designation of span number 
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Fig. 2-2 End point designation 
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V. SOLVING A SIMULTANEOUS EQUATION USING COMPUTOR 709 O 



Before solving a simultaneous equation ve decide to solve two 



different loading conditions: (l) favorable loading condition (see Fig. 



3 )^ ( 2 ) most critical loading condition (see Fig. . 

For the girder deflection ve must also think about shear defor- 
mation because of having relatively short span length and deep depth. 
For this purpose ve can find deflection equations of simple beam at the 
condition of both ends fixed and concentrated at some span point. 

One can use the next simple equation for the concentrated load 
at the center of the beam. 



W 





When the concentrated load is acting at the point of l/k span 



length; deflection can be expressed as follows: 



W. 



P 5 




Where W^. . .Deflection at the point of load applied 
A-j^. . .Web area of b beam 



P . . .Concentrated load 



For the computed deflection; see curves 4-1; 4-2 and 4-3- 
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FIG. 4-1 UNIT LOAD DEFIECTION OF LONGITUDINAL CENTER GIRDER 
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FIG.4-2UNIT LOAD DEFLECTION OF DECK CENTER GIRDER 
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FIG.A*3 unit load deflection of longitudinal BHD 
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One of the computed results is tabulated below as Table 5 * 



X 


Deflection angle in Radian 
X 10-4 


X(l) 


0 . 475742 


X( 2 ) 


-0.413502 


x(3) 


-0.226263 


x(4) 


0.128071 


x(5) 


0.171848 


x(6) 


-0.369586 


x(7) 


-0.501310 


x(8) 


0.550122 


x(9) 


0.647280 


x(io) 


-0.336503 


x(ii) 


0.028208 


X( 12 ) 


-0.080296 


X( 13 ) 


-0.070333 


x(i4) 


0.096957 



Table 5 • Acquired unknown variables X at the 
forward frame (load l) in two 
dimensional analysis. 

From these values we can get moment at edge 1 and 2 (left and 
right side or upper and lower side) for every span. You can see one of 
the results in Table 6 . 
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SPM No. 


MOMENT 1 X 103 


MOMENT 2 X 103 


1 


-0.242163 


0.914909 


2 


-0.894309 


1.194515 


3 


-1.177446 


0.960957 


k 


-0.933445 


10.523114 


5 


0.388150 


-0.975332 


6 


0.936319 


-1.158020 


7 


1.201756 


-0.847029 


8 


0.799642 


-11.777653 


9 


-0.059339 


0.242128 


10 


0.039037 


-0.020693 


11 


-0.043852 


-0.170775 


12 


0.047344 


-0.027552 


13 


12.349345 


-10.523304 


14 


-0.169027 


-0.328803 


15 


-0.283449 


-0.571738 



Table 6. Moment of the forward frame 

From these end moments we can get very easily a middle portion 
moment by simply superimposing loading to each span. This work can be 
done very easily if we use the Computer for every span. The results 
are plotted in figures . 

VI. FIXED POIET METHOD 
A. General 

For the purpose of checking above results and investigating 
shear effect^ we will apply here numerical computation by the method of 
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fixed point. The fixed point method is originated from the three moment 
equation. In this approach^ we have the following assumptions: 

(1) All nodal points are fixed in space. 

(2) No shear force is considered. 

We will not discuss here in detail about this derivation (see Appendix 
but will apply this method in practical calculation using the same load- 
ing conditions to the same scantling frame ship. 
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Then the moment at the ends of the span can be expressed by the follow- 
ing formulas . 



M 



ef 



M 



de 




I *2 



4-2 



- R 
e e 



A.. 



R. 



?0e 
"e 7 ^ 
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Here 

/• 0^ p 



is the span length 
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In generally fixed point 

ze - (5n 



^n = 



L 



6e ^ + 3 - 6^) 



4 = 



Y.(x) -[})y 






L 



6e -^ + 3 (2w - w„) 



and stiffness 



«n - -r- 

On 



1 

I ^ 7 — 

n <^n 



2 - 3 



A^y 



similarly 



A 



^n = “2^ 



1 - 



2 - 3 



^n 



And Rg and can be found out in some handbooks » 

That is: 

L e = AqX 
= A^d-x) 

Aq is the area below the simple moment curve, 
and X is the distance from the left side to the center 
of the area Aq . 

For our computation purpose we can find out the following data from the 
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In generally fixed point 

ze - t-n 



K = 



6E-^+ 3 - 6^) 



A = 



Ti (xJ ^ CaJ]q 



tr 



6e -^ + 3 (Ew - (On) 

^-v 



n 



and stiffness 






1 - 



n 2 _ 3 






similarly 



A 






1 - 



2 - 3 






And Rg and Lg can be found out in some handbooks. 

That is: 

Le = AgX 

R^= Ag(l-x) 

Aq is the area below the simple moment curve, 
and X is the distance from the left side to the center 
of the area Aq . 

For our computation purpose we can find out the following data from the 




handbook. 
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If once we get moment for each span^ then we can compute each joint 
moment using stiffness and fixed points by one-by-one adding process as 
follows : 




But in doing so^ we must be very careful in the sign for some opening 
span. (Ref. illustrative calculation of moment distribution Table 12.) 



Now here is the practically computed results listed. 







^in 


^n 

^n 


^15 

Unit X -7— 1 - 

^15 


1 


30828.851 


280.539 


214.955 


2.480 


2 


847 ij-l . 710 


301 . 006 


281.528 


3.248 


3 


lj.5i97.227 


137.987 


327.546 


3.779 


k 


6081j-7.290 


137.987 


440.963 


5 .088 


5 


5282 . 347 


301.006 


17.548 


0.202 


6 


25839.827 


125.984 


205.104 


2.366 


7 


21014.376 


125.984 


166 .801 


1.924 


8 


3656.705 


301.006 


12.148 


o.l 4 o 


9 


25839.827 


125.984 


205.104 


2.366 


10 


20716.733 


125.984 


164.439 


1.897 


11 


3083.156 


301.006 


10.242 


0.118 


12 


30828.851 


198.454 


155.345 


1.792 


13 


23604 .205 


198.454 


118.940 


1.372 


Ik 


22181.710 


301 . 006 


73.691 


0.850 


15 


24303.928 


280.539 


86.632 


1.000 



Table 7 
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Define K = 6E ^ 

^15 

In the next table the computed values of , ^n/^ ? 

and Wj^/k are shown • 





\i/ ^n 


^n/K 


^ /C 
n/ 


"^n/K ' 


1 


0.232 


1.055 


0.236 


1.022 


2 


0.219 


1-375 


0.224 


1.382 


3 


( 0 . 250 ) 0.248 


1.419 


0 .256 


1.428 


h 


0.200 


2.159 


0.203 


2.164 


5 


0.240 


0.813 


0.245 


0.818 


6 


0.223 


1.107 


0.221 


1.104 


7 


(0.250)0.275 


1.154 


0.277 


1.189 


8 


0.240 


(3.816)3.021 


0.205 


2.920 


9 


0.327 


0.561 


0.262 


0.516 


10 


0.327 


0.077 


0.328 


0.077 


11 


0.327 


0.092 


0.328 


0.092 


12 


0.327 


0.135 


0.328 


0.133 


13 


0.260 


1.969 


0.222 


1.893 


Ik 


0.266 


0.610 


0.268 


0.612 


15 


0.284 


0.703 


0.285 


0.703 



Table 8 



On the next page, 



f ! P' and A/ a' tables are also prepared. 
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K/K 


1 


0.308 


0.302 


2 


0.288 


0.208 


3 


0.341+ 


0.329 


k 


0.254 


0.250 


5 


0.324 


0.315 


6 


0.283 


0.287 


7 


0.383 


0.379 


8 


0.257 


0.315 


9 


0.355 


0.485 


10 


0.424 


0.485 


11 


0.424 


0.485 


12 


0.488 


0.485 


13 


0.285 


0.351 


li+ 


0.366 


0.362 


15 


0.398 


0.396 



Table 9~1 
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The computed end moment by its own span load are listed in 
Table 10 . (We restricted our computation in the foiTward frame only.) 



SPAN No. 


LOAD 


I 


LOAD 


II 


MOMENT 1 X 103 


MOMENT 2 X lO^ 


MOMENT 1 X 103 


MOMENT 2 X 103 


1 


0.42470 


0.52657 


0 . 42470 


0.52657 


2 


0 . 47568 


0.48144 


0.47568 


0 . 48144 


3 


0 . 81965 


0.79450 


0.81965 


0.79450 


k 


1.03658 


1.03886 


1.03658 


1.03886 


5 


-0.44029 


-0.54524 


0 


0 


6 


-0.46159 


-0.49807 


0 


0 


7 


-0.87103 


-0.89650 


0 


0 


8 


-0.95682 


-0.13608 


0 


0 


9 


0 


0 


0 


0 


10 


0 


0 


0 


0 


11 


0 


0 


0 


0 


12 


0 


0 


0 


0 


13 


-7.95302 


-5.59423 


-7.95302 


-5.59423 


14 


0 


0 


0 


0 


15 


0 


0 


-6.82683 


-6.89544 



Table 10 



From the above computed end moment at edge 1 and 2 we can dis- 
tribute moment by using stiffness and fixed points and , (See 
Table 12 .) Then we can get the following final end moment. 
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Table 11. Final end moment in tons-inch 





LOAD 


I 


LOAD 


11 


SPAN No. 












MOMENT 1 X 103 


MOMENT 2 X lo 3 


MOMENT 1 X 103 


MOMENT 2 X lo 3 


1 


0.2968 


0 . 7627 


0.3396 


0.7616 


2 


0.7425 


1.1749 


0.7536 


1.1315 


3 


1.2091 


0 . 4132 


1.1643 


0.5683 


4 


0 . 3464 


6.2266 


0.5379 


9.6970 


5 


-0.3354 


-0.7174 


0.1312 


0.2012 


6 


-0 .6946 


- 1-3596 


0.2248 


-0.8502 


7 


-1.4136 


-0.1739 


-0.9191 


2 . 4365 


8 


-0.0980 


-7 • 3683 


2.6901 


-1.9521 


9 


0.1419 


0.2698 


-0.0935 


0.3396 


10 


-0.0228 


-0.0202 


-0.0236 


-0.0080 


11 


0.0535 


0 . 0342 


0 . 0689 


-0.0328 


12 


-0.0759 


-0.0668 


-0.2536 


-0.0304 


13 


- 8.4853 


-6.2266 


-10.2938 


-5.6970 


14 


0.0662 


-0.1935 


0.0049 


0.0377 


15 


0.6203 


-1.1170 


-6 . 3110 


-8.3417 
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Table 12 . Illustrative calculation of moment distribution 



for span No. 1 


load 


-25.3 63.6 


150.7 - 424.7 


-87.1 




27.4 


-12.3 


-5.1 


0.6 


22.3 


-11.7 26.5 


26.8 




7.2 




34.0 


3.8 


-7.7 


- 0.6 


1.9 


-0.1 


-1.3 


3.1 -9.0 


-7.1 






-7.1 




-1.9 


2.7 


-1.2 


- 0.4 


-10.2 


2.3 


-1.2 




0.2 2.5 


0.8 


-1.0 


2.5 




3.3 




2.5 




-0.8 




1.7 




-0.3 , ^ 


0 c- -10.1 


1.0 


3.5 


0.3 


-0.5 


1.3 


3.0 



+ 




(c) 




(a) 



- 



J 



+ 



-h 



424.7 



526.5 

500.0 



1 4 0.0 

131.0 



-43.0 

-40.5 



10.1 
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Once ve get each span*s end moment then ve can get various 
point moments very easily by superimposing process. And these numeri- 
cal computational results are plotted in figures . 

Also stress curves are available in this present paper^ but we 
did not compute the moments and stresses for all the frames by the 
fixed point method. We only calculated the forward frame for both 
loading conditions . 

The reason is this: the fixed point method is^ of course^ a 

two dimensional one and each nodal point is fixed in space; therefore^ 
a little difference of water head will bring us a very insignificant 
computational difference. 

VII. DISCUSSION 

Most of the naval architects are anxious to learn about the 
three-dimensional effects on the transverse frame strength of the ship 
But the computation is rather lengthy and tedious to work out so the 
architects estimate the effect on transverse frame strength by the 
simplified method which satisfies them. In this paper we dealt with 
this lengthy calculation by use of the computer. It cannot be done 
without this machine because so many unknown variables are interacted. 

Even a casual glance at the graph will show that there are not 
many differences between the two-dimensional and the three-dimensional 
analysis . 

Especially in load condition |L; the results are nearly the same 
Also in load condition 2^ we do not see much difference between these 
two . 

In load condition 1^ at the forward frame (Fig. 5) can see 
also almost no differences between the fixed point method and our 
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present method in the upper half of the frame that is span No. 2; 5^ 

6; 9; 10^ 11 and l4. That means that lightly loaded spans do not have 
any shear effect at all^ hut loading is steadily increased^ then the 
discrepancy is noticeable. 

We are now going to examine these graphs individually for more 

detail . 

A. Load Condition 1 

(1) This loading condition is ballasted with sea water in 
the center tank and the side shell is exhibited to the hydrostatic load. 
This is one of the most favorable loading conditions. In this loading 
condition^ the forward frame does not show us much discrepancy between 
two-dimensional and three-dimensional analysis. In this loading condi- 
tion^ junction points 1-2 and 2-1; 2-2 and 3"1^ 3“2 and 4-1; also 

and 6-1; 6-2 and T~l; and and 8-1 remain in space because of the 
symetric loading condition. Therefore the fixed point method and our 
slope-deflection method can be comparable in quantities . 

ThuS; this means if we neglect shear effect, we will lose al- 
most 4,000 Tons-inch moment at the span edge which is most heavily 
loaded. 

(2) Middle frame: We can see a little difference between 

these two methods, but not too much. This means that the middle frame 
is the most effected frame among the three frames because the girders 
will deflect most at the center. 

(3) In the aft frame, we can see almost no effect of the 
three-dimensional effect on the moment curve. Almost the same as the 



forward frame . 



34 . 



B. Load Condition 2 

(1) We think this loading condition is most critical^ but 
still two-dimensional analysis and three-dimensional analysis have the 
same aspect. They don*t change very much; but here are several inter- 
esting aspects being displayed. That iS; in our method we fixed joint 
1-1 and 9~2 and 4-2 and 13-2 so that as a whole side frame and BHD bent 
something like a semi -fixed beam at the bottom and deck; but by the 
fixed point method nodal points are fixed in space so that they have 
almost the same aspect as in load condition 1. But from this curve we 
can see also a big difference in heavily loaded lower spans . 

(2) In the middle frame we can see some conceivable dis- 
crepancies between those two methods; but still in good shape; in 
three-dimensional analysis; however; the moment at point 15-1 and 9-2 
is reversed in direction because of the restraints by girders released 
at 15-1 and 13-1 points. The other spans have almost the same results; 
nevertheless one can see a little higher moment in the frames 9; 10; H; 
12 and l4. 

(3) In the aft frame; the aspect is almost the same as in 
the forward frame. We cannot see any noticeable discrepancies between 
those two approaches. 

CONCLUSION 

Most designers of super tankers are worried about the three- 
dimensional effect on the transverse frame strength. But from our com- 
putation we cannot find large discrepancies between two-dimensional and 
three-dimensional analysis. However; in the three-dimensional analysis; 
calculation is more detailed and tedious than in the two-dimensional. 
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For our present discussion; three-dimensional computation is 
more academic than practical; othervise we would not use any simplified 
design process such as Professor Schade*s Design Curves for Cross 
Stiffened Plating. 

Thus our conclusion is: forget about the three-dimensional 

analysis; unless you intend to use any simplified design process pre- 
pared for the three-dimensional analysis; but take into account the 
shear effect on the two-dimensional analysis on super tanker designs. 
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FIG. 5 FORWARD FRAME MOMENT CURVE (LOAD 1) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 

JIXED POINT METHOD 




• .*1 
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FIG. 6 MIDDLE FRAME MOl^NT CURVE (LOAD 1) 

T\70 DIMENSIONAL 

THREE DIMENSIONAL 
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FIG. 7 AFT FRA14E MOMEOT CURVE (LOAD i) 

TWO DIl-lENSIONAL 

THREE DIMENSIONAL 
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FIG. 8 FORWARD FRAME MOl^lENT CURVE (LOAD 2) 
T\^0 DIMENSIONAL 



THREE DIMENSIONAL 

FIXED POINT METHOD 
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FIG. 9 MIDDLE FRAME MOMENT CURVE (LOAD 2) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 




2000 IN-TONS 




I 




FIG. 10 AFT FRAME MOMENT CURVE (LOAD 2) 



TTTO DIMENSIONAL 

THREE DIMENSIONAL 




2000 IN-TONS 
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FIG. 11 FORWARD FRAME STRESS CURVE (LOAD 1) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 

FIXED POINT METHOD 
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FIG. 12 MIDDLE FRAME STRESS CURVE (LOAD 1) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 
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FIG. 13 AFT FRAME STRESS CURVE (LOAD 1) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 
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FIG* 14 FORWARD FRAME STRESS CURVE (LOAD 2) 



DIMENSIONAL 

THREE DIMENSIONAL 

FIXED POINT METHOD 
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FIG. 15 MIDDLE FRAME STRESS CURVE (LOAD 2) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 




\ 

\ 

\ 

\ 

\ 



47 . 



FIG. 16 AFT FRAME STRESS CURVE (LOAD 2) 

TWO DIMENSIONAL 

THREE DIMENSIONAL 
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IX . APPEM)IX 



1. The three moment equation 



or 





Then by using the second moment area theorem to obtain 



tLc and 






Mxdx 



Therefore simply 



^ 3 \ " I Vc I 

■ ■ ^ I ^”a 3 ^ + I ^“c I '-r) 



Simplifying the above expression, the three moment equation is: 

b.«L " " A " ^ «C " 5 




6Ar 
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2. Fixed point method 




By three moment theory (Ref. Appendix l) first three support (from left) 






- 0 



Then 



(where = 0 ) 

set i , etc. 



0 + 2M-Kp (^b + ^c) + Med = 0 



Ml 



^bc _ 






“od "«b " «c) 



By geometry 



Therefore finally 



^bc A^c 



A^c 



2(f^ + 



Apply three moment equation to span C and D 

Ac«c b 2M^d «c + (-if * Hie » 



using relationship 



Mbc = - ^ 



M. 



cd 



A., 



Med (2^^ + 2e^ - e^) - M^e ^d= 0 



^d ^ £d 

A-d' p/? + A/c 

^ Oc 






n 



In general 




•n - 1 



NOTE: The position of zero moment is a function of geometry and the 
degree of fixity at the ends and independent of the type of loads or 
the magnitude of loads . 



Apply three moment equation to span D and E is unloaded, E is 
loaded) . 



Where = Moment of field load moment area about right 
support of span E 

= Moment of field load moment area about left 
support of span E 



LOADED SPAN E 




A 








Eliminate M^^ and M^^ by using: 



51 - 



Then we obtain 

(2 



A., 



^ U * M.^c + 6 ^ e„ . 0 



de '• d e 






ef e 



2 ''e 



«ef (2^e + 2^f - 7^ ^) + «ae4 + « 7 ^ 4 ' 



= 0 



from general equation 

A. 



Ai 



2^d + 2^ - ^d - '-e 



F 






e 



2^g + 2(f|. - t'f - 



Then above equation becomes 

A' 



Mde 'i: - ^ 6 



e.2 



p' 

> o 



M. 



bf + 6 ^ 



= 0 



= 0 



Finally we get 



M. 






\ 



e ^'e 



M 



de 



f 

\ ^e ^ 



^2 

6 

^p>2 



L -R ^ 
e e /Ve 

R L ^ 

^e -^e />f 



3 . Stiffness factor 

Stiffness factor is defined as the reciprocal of the angle of 
rotation caused by a unit moment. 

From the elementary mechanics of materials we know the 
following relation between angle of rotation and moments 




A = ^ L 
El 6 






M 



B L 



El 3 
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By these relations we can easily obtain 

t ^ 

" 3El' ^CD 



M = 1 , 
CD 






a. 



M. 



a; CD 



taxi<^ 



t. 



Aq 



CD 3EI( 



A'c 6 ei^ 



If we use the notation t ^ to denote the stiffness of the left 
hand side member, and substitute by ( ) we will get 



tanizi 26 cA’c-Ac^c 2 ^c( 



CD 



6EI^( ^c-^c) 



6ei,(6c-Ac) 



6ei, 



^c-^c 






In general; the stiffness of the left hand member at the right 
hand support t 



= 6e hi 



1- An 

^ 

■^n 2 - SAn 






Similarly we can express in terms of right hand member 

-t&ni = ^d , - 'hd -Ad)^(i 

CD 3EI(3^ 6EI^ 6EI(J^A^^ 

_ ^d '2-td. - -3^d + 3 Ad 
6Eia 



Therefore 



ta 



tc 6ei^ 
Ka = 






:^ + 3 






6e 






6 
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In general 




members* stiffness^ then we have 






6e -^ + 3(2^ - t^) 

*'n 

Similarly we can apply the above formula from right side to 
left side. The results are: 



= 6e 






1 



n 



'n 2 _ n 



A = 



W- OJ 



n 



6e " + 3(2w - 






n 



■21 



no o 

^ ^'tr. 



'^n 



Zt - t 



n 



I 



6e -^ + 3(z^ - t ) 



5 ^. 
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